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Abstract 

We demonstrate how measuring orbital period deriva- 
tives can lead to more accurate distance estimates 
and transverse velocities for some nearby binary pul- 
sars. In many cases this method will estimate dis- 
tances more accurately than is possible by annual 
parallax, as the relative error decreases as Un- 
fortunately, distance uncertainties limit the degree 
to which nearby relativistic binary pulsars can be 
used for testing the general relativistic prediction of 
orbital period decay to a few percent. Neverthe- 
less, the measured orbital period derivative of PSR 
B1534-I-12 agrees within the observational uncertain- 
ties with that predicted by general relativity if the 
proper-motion contribution is accounted for. 
Subject headings: Gravitation - Relativity - Pulsars; Gen- 
eral 

1. Introduction 

The most common method for determining the dis- 
tances to radio pulsars is based on their dispersion 
measure and models of the Galactic d istribution of 
free electrons (Taylor & Cordes 1993). These dis- 



tance estimates typically have an uncertainty of 30%. 
Distances may also be determined by measuring an- 



nual parallax, based on either timing (Ryba & Taylor 



1991) or interferometric measurements (Gwinn et al 



1986). Hi absorption by interstellar hydrogen is also a 
common distance indicator (Frail & Weisberg 1990). 
However, no pulsar has a distance estimate more ac- 
curate than ^ 10%, and for all but two, the errors are 
greater than 20%. 

Any acceleration of a pulsar along the line of 
sight will change the observed pulse period deriva- 



tive P. As Shklovskii (1970) pointed out, an ap- 
parent acceleration occurs when the proper motion 
is significant. The magnitude of this contribution is 
P(pm)/P = v^/{cd), where P is the pulse period, v is 
the transverse velocity, d is the pulsar distance, and 
c is the speed of light. For many millisecond pulsars, 
this effect is of similar magnitude to the intrinsic pulse 
period derivative, making it hard to determine accu- 
rately from timing data either the intrinsic pulse pe- 



riod derivative (Camilo, Thorsctt & Kulkarni 1994) 



or the distance and transverse velocity. 

2. Distances and Velocities 

This apparent acceleration also applies to orbital 
period derivatives, and the contribution is Pb('pm) I Pb = 



V /(cd), where Pb is the orbital period (Damour & 



Taylor 1991). In fact, for many nearby millisecond 
pulsars, it is expected to completely dominate future 
observed orbital period derivatives. This means that 
v'^/d can be obtained and, when combined with the 
measured proper motion /i = v/d, the distance and 
transverse velocity can be easily separated. Hence, 
the proper-motion contribution to the pulse period 
derivative can also be determined, giving accurate es- 
timates of the intrinsic pulse period derivative and 
hence the magnetic field strengths, ages and spin- 
down luminosities of binary millisecond pulsars. 

The amplitude and functional form of the residu- 
als from a least-squares fit to the observed pulse ar- 
rival times, if one parameter is set to zero, is often 
called the "timing signal" for that parameter. For 
proper motion, the timing signal is often relatively 
large, with its amplitude increasing linearly with time. 
With continued measurement therefore, its relative 
error decreases as t^^-^. The peak-to-peak amplitude 
AT(p„j) of the timing signal owing to the contribution 
of the proper motion to the observed orbital period 
derivative, is 



a sin i 2tt v 
c Ph cd 



(1) 



where a is the semi-major axis of the pulsar's orbit 
and i is the orbital inclination. The accuracy of dis- 
tances obtained in this way are limited by the accu- 
racy of the orbital period derivative measurements. 
Their accuracy, and therefore the accuracy of dis- 
tances improve as t~'^^^ . The fact that the relative 
error in both of these critical parameters decreases 
in such a spectacular fashion with time demonstrates 
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the power of this method for determining distance and 
transverse velocity. 

Table |l| shows the predicted size of the timing sig- 
nal ATpm after 10 years of regular timing observa- 
tions for a selection of binary millisecond pulsars. 
Where the proper motion was not available, we used 
the median transverse velocity for millisecond pul- 
sars of 69 km s~^ was used. Also shown is the tim- 



ing signal due to parallax (Ryba & Taylor 1991), 
AT^ = cos^ 6 / {2cd) , where r is the radius of the 
Earth's orbit and 6 is the angle between the line of 
sight to the pulsar and the Earth's orbital plane. Af- 
ter 10 years, the new method will provide better dis- 
tance estimates than parallax measurements. This is 
possible because ATp^ oc , while AT^ is constant. 
If an rms timing residual of 1.0 /iS could be obtained 
it would be possible to determine distances this way 
for several of the currently known binary millisecond 
pulsars. On average, for the pulsars listed in Table 0, 
2-3 years of precise timing data have been recorded 
by various observers. So, to reap the rewards of this 
method a further 7-8 years of precise timing will be 
required. 

3. Other Contributions to Period derivatives 

Many other effects could contribute to an observed 
orbital period derivative; for example, changes in the 
gravita tional constant jDamour fc Taylor 1991 ), tidal 
effect s ( Will 199^, Arzoumanian, Fruchter fc Taylor 



1994, Applegate fc Shaham 1994), companion mass 



loss ( [ben &: Tutukov 1986) and accelerations in glob- 



ular cluste r potentials ( Blandford, Romani fc Apple- 
gate 1987 ) . These contributions are indistinguishable 
from the proper-motion contribution, and so it is im- 
portant to determine which of them are significant. 
Known pulsars possess one of 5 types of compan- 
ion: a neutron star, a main sequence star, a white 
dwarf, a very low mass helium star, or a planetary 
system (Wolszczan 1994). Fortunately, binary pul- 
sars with either white dwarf or neutron star compan- 
ions are very "clean", and their orbital periods not 
affected by tidal or mass-loss effects ( p^ben fc Tutukov 
1986). Systems with low-mass companions such as 



PSR B1957-f20 possess large orbital period deriva- 



tives, possibly caused by tidal effects (Applegate & 
shaham 1994). The small asini induced by planetary 
companions in the pulsar orbit makes it extremely 
difficult to measure their orbital period derivatives. 
The only significant contributions to the orbital pe- 



riod derivatives in neutron star and white dwarf sys- 
tems are those due to acceleration in the Galactic 
potential Pf,(fc2) , Galactic differential rotation P^dr) j 
proper motion Pb{pm)i a-nd general relativity P^j^gj-y 
Table |l| lists those contributions showing, that the 
proper-motion term will dominate for many of the bi- 
nary millisecond pulsars. 

For the nearby millisecond pulsar J0437— 4715, the 
uncertainty in Pi)(kz) is approximately 1% of Pb(pm)- 
Hence, measurement of the orbital period derivative 
will ultimately provide a distance estimate which is 
limited in accuracy to about 1%. If the distance could 
be independently estimated with superior accuracy, 
it would be possible to determine the acceleration of 
the binary in the Galactic gravitational potential and 
thereby constrain the distribution and composition 
of dark matter ( Flynn fc Fuchs 1994 ). However, pul- 
sars such as PSR J2317-f 1439 with large z-heights 
are probably better suited to such an exercise, be- 
cause contribution from the Galactic acceleration in 
such pulsars is comparable to the contribution from 
the proper motion. This emphasises the importance 
of monitoring known binary millisecond pulsars and 
searching for new ones. 

4. Implications for Tests of General Relativ- 
ity 

The double neutron-star system PSR B1534-I-12 
has been predicted to provide an even better rela- 
t ivistic laboratory t han the binary pulsar B1913-f 16 
( Arzoumanian 1995 ). Unfortunately, the distance to 
this pulsar is known only to an accuracy of some 30% 
( Taylor fc Cordes 1993| ), and therefore the proper 
motion contribution to the orbital period derivative 
is uncertain by a similar amount. Recent measure- 
ments (Arzoumanian 1995) indicate that the pre- 
dicted orbital period derivative due to gravitational 
wave emission Pb(gr) is —1.924 x 10"^'^, whereas the 



observed value P; 



b{obs) 



is only -1.5 ± 0.3 x 10 



-13 



Using the dispersion-measure distance of 0.68 ± 0.2 
pc and the measured proper motion, the contribu- 
tion to the observed value from the proper motion 



is Pbiprn) = 0.40 ± 0.12 X 10-". 



Since Pi 



P> 



b{obs) 



P> 



b{pra) 



b(gr) — 

-1.9 ± 0.3 X 10"", the observed 
value is in excellent agreement with the general rela- 
tivistic prediction. Unless the distance estimate can 
be improved, the orbital period decay due to the emis- 
sion of gravitational waves cannot be verified to bet- 
ter than ~ 5% in the PSR 1534+12 system. This is 
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a surprising rc^siilt. wliic;h underlines the importance 
of obtaining independent distance estimates to this 
system. 
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Pulsar 
Name 


AT 


AT, 


Contributions to Pb/Pft(ohs) (xlO 

Pb{kz)/Pb Pb{dr)/Pb Pb{gr) / Pb 


-19 g-l) 
Pb{pm)/Pb 


J1713+0747^ 


0.3 


1.0 


-0.89 


-0.29 


- 


0.88 


BI855+O92 


0.5 


0.8 


-0.029 


-0.002 


-0.001 


0.96 


J0613-02003 


1.5 


0.4 


-0.33 


-0.67 


-0.25 


2.4 


J2317+1439^ 


1.8 


0.5 


-2.2 


0.89 


-0.043 


2.8 


J2145-07505 


2.2 


2.4 


-1.3 


0.033 


-0.009 


2.1 


J0751+1800^ 


2.8 


0.6 


-0.98 


-0.77 


-21.8 


2.6 


JOO34-O5345 


3.4 


1.2 


-2.8 


0.17 


-0.15 


5.3 


B1913+16^ 


3.8 


0.1 


-0.064 


5.4 


-845.3 


0.74 


J2019+2425^ 


4.5 


0.7 


-0.11 


0.45 


- 


12.6 


JIOI2+53O79 


7.0 


1.5 


-1.7 


0.16 


-1.8 


10.1 


J1022+021'' 


13.5 


2.0 


-1.8 


0.048 


-0.010 


8.7 


J1455-33303 


22.4 


1.5 


-0.69 


-0.24 


- 


73.1 


J0437-4715" 


28.5 


1.2 


-0.57 


0.058 


-0.004 


67.5 


B0655+642 


31.5 


1.4 


-0.65 


-0.16 


-2.1 


10.9 


BI534+1212 


70.5 


1.3 


-1.8 


-0.24 


-52.8 


11.0 



Table 1: Predicted orbital period derivatives and timing signals. ^Camilo, Foster & Wolszczan 1994 
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